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Abstract 

To investigate the impacts of biophysical factors on light response of net ecosystem exchange (NEE), C0 2 flux was measured 
using the eddy covariance technique in a winter wheat field in the North China Plain from 2003 to 2006. A rectangular 
hyperbolic function was used to describe NEE light response. Maximum photosynthetic capacity (P max ) was 46.6±4.0 umol 
C0 2 m~ 2 s -1 and initial light use efficiency (a) 0.059±0.006 |imol ixmol -1 in April — May, two or three times as high as those 
in March. Stepwise multiple linear regressions showed that P max increased with the increase in leaf area index (LAI), canopy 
conductance (g c ) and air temperature (7" a ) but declined with increasing vapor pressure deficit (VPD) (P<0.001). The factors 
influencing P max were sorted as LAI, g c , T a and VPD. a was proportional to In(LAI), g c , T a and VPD (P<0.001). The effects of 
LAI, g c and T a on a were larger than that of VPD. When I a >25°C or VPD>1.1 -1.3 kPa, NEE residual increased with the 
increase in 7" a and VPD (P<0.001), indicating that temperature and water stress occurred. When g c was more than 14 mm 
s -1 in March and May and 26 mm s -1 in April, the NEE residuals decline disappeared, or even turned into an increase in g c 
(P<0.01), implying shifts from stomatal limitation to non-stomatal limitation on NEE. Although the differences between 
sunny and cloudy sky conditions were unremarkable for light response parameters, simulated net C0 2 uptake under the 
same radiation intensity averaged 18% higher in cloudy days than in sunny days during the year 2003—2006. It is necessary 
to include these effects in relevant carbon cycle models to improve our estimation of carbon balance at regional and global 
scales. 
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Introduction 

Vegetation productivity is the foundation of carbon sequestra- 
tion and grain yield formation. For the limit on field observation 
techniques, early studies mainly focused on photosynthesis at the 
leaf level. The biochemistry and ecophysiology of leaf photosyn- 
thesis have been well understood and parameterized [1,2]. 
However, it is difficult to measure leaf photosynthesis over long 
periods and upscale photosynthetic rate from the leaf level to the 
canopy or ecosystem level because of the non-linear distribution of 
leaf area and radiation intensity within the vegetation canopy. 
With the development of micrometeorological techniques, espe- 
cially the eddy covariance (EC) method, net photosynthetic rate of 
vegetation could be direcdy measured at the ecosystem level. The 
EC technique has been widely used in C0 2 flux measurements in 
forest, grassland and farmland ecosystems [3-9] . 

Plant photosynthesis is primarily driven by incident solar 
radiation. Light response models, including rectangular hyperbolic 
model [10,11] and non-rectangular hyperbolic model [12,13], 
were developed to describe the relationship between daytime net 
ecosystem C0 2 exchange (NEE) and photosynthetically active 



radiation (PAR). Besides solar radiation, factors influencing 
daytime NEE include environmental variables such as air 
temperature, vapor pressure deficit and soil water content and 
biological variables such as leaf area index and canopy conduc- 
tance. These factors may be considered by (1) estimating the bias 
of simulated NEE as functions of influencing factors [14,15], or (2) 
revising light response parameters as functions of influencing 
factors [9,11,16-19]. The effects of biophysical factors on light 
response parameters are often estimated using stepwise multiple 
linear regression models [9] . However, these models have not been 
used to assess the influence of variables on NEE residual so far. 

Long-term observations and simulations have shown a world- 
wide decrease in surface solar radiation (global dimming) from 
1950s to 1980s, with a partial recovery (brightening) in 1990s in 
some areas (e.g. high and middle latitude of the Northern 
Hemisphere) [20]. During the dimming period, direct radiation 
declined remarkably, whereas diffuse radiation enhanced. The 
increase in diffuse fraction of surface solar radiation may be 
attributed to increasing aerosol and/or cloudiness, as both factors 
tend to enhance scattering in the atmosphere [20] . Light responses 
of photosynthesis under various sky conditions were investigated in 
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forests [8,21-23] and croplands [24,25]. For the tall and dense 
vegetations, photosynthesis may be enhanced by diffuse radiation 
because diffuse light distributes more effectively within the canopy 
compared with direct light [26]. After considering the effects of 
direct and diffuse radiation on canopy photosynthesis in the 
model, Mercado et al. [27] found that the increases in diffuse 
fraction enhanced the global land carbon sink by 23.7% during 
the period from 1960 to 1999. 

China is the largest wheat producer and consumer in the world 
[28] and continuously attempts to increase its production to ensure 
national food security. As one of large food production regions in 
China, the North China Plain (NCP) produces about half of the 
country's wheat [29]. In this study, C0 2 flux was measured 
continuously using the EC technique in a winter wheat field in the 
North China Plain for 4 years. The objectives are to (1) investigate 
NEE light response and the influencing factors, and (2) assess the 
effects of sky conditions on NEE light response. This study will 
improve our knowledge on the parameterization of carbon cycle 
models and the scenario analyses of carbon sink in the future 
under changing climate. 

Materials and Methods 

Study site 

This study was conducted at Yucheng Comprehensive Exper- 
iment Station, Chinese Academy of Sciences (36°57'N, 1 16°38'E, 
23.4 m). It is located at the North China Plain, with a temperate 
monsoon climate. Mean annual temperature is 13.1°C and annual 
solar radiation is 5242 MJ m . Annual precipitation is about 
528 mm. Soil organic content is 1.21% and pH value is about 7.9. 
The typical cropping system in this region is the biannual rotation 
with winter wheat and summer maize. In this study, winter wheat 
was planted in mid/late October and harvested in early/ mid June. 
The detailed field management was described by Tong et al. [30] . 

Field observations 

C0 2 and latent heat fluxes were measured by the eddy 
covariance system with a three-dimensional sonic anemometer 
(model CSAT3, Campbell Sci. Inc., USA) and an infrared open- 
path C0 2 /H 2 0 gas analyzer (model LI-7500, Li-Cor Inc., USA). 
The eddy covariance system, mounted at the height of 2. 1 m, was 
used to measure 3-D wind speed, air temperature, humidity and 
C0 2 concentration above the canopy. Raw data were collected at 
10 Hz and recorded by a CR5000 datalogger (model CR5000, 
Campbell Scientific Inc., USA). The CR5000 datalogger can store 
data and programs on a PC card. The card can be carried to a 
computer and the computer reads PC cards via its PCMCIA card 
slot direcdy. 

Anemometers (model A100R, Vector, UK) and psychrometers 
(model HMP45C, Vaisala, Finland) were installed at heights of 
2.05 and 3.25 m above the ground. Photosynthetically active 
radiation (PAR) was measured using a quantum sensor (model Li- 
190SB, Li-Cor Inc, USA). Solar radiation and net radiation (ffj 
were measure by a pyranometer (model CM11, Kipp & Zonen, 
Delft, The Netherlands) and a net radiometer (model CNR-1, 
Kipp & Zonen, Delft, The Netherlands), respectively. Two soil 
heat flux plates were buried in the depth of 2 cm, one between 
rows and another between plants. Soil temperature was measured 
at 0, 5, 10, 30 and 50 cm depths. Soil water content (SWC) at the 
depths of 20 cm and 30 cm was measured with time domain 
reflectometers (TDR) (model CS616, Campbell Sci. Inc., USA). 
Rainfall was measured with a rain gauge (model 52203, Rm 
Young MI, USA). All meteorological data were recorded with a 



data logger (model CR23x, Campbell Sci. Inc., USA) and were 
stored at intervals of 30 min. 

Biomass, leaf area index (LAI) and plant height were measured 
every 5 days during the growing season of winter wheat. There 
were three sampling plots (replications) for each measurement. 
Twenty plants were sampled continuously for each plot and the 
plots were selected randomly. Leaf area was measured with Leaf 
area instrument (model Li-3100, Campbell Sci. Inc., USA). In this 
study, the main growing season began when the cropland turned 
from the carbon source to the sink (5-day moving mean C0 2 flux 
from positive to negative), and it ended when shifted reversely. 

Flux data quality control 

C0 2 and latent heat fluxes were calculated as follows: 



F c = p(w'c') 



XE = Xp(w'q') 



(1) 



(2) 



where p is air density, w' the vertical wind velocity, c' C0 2 
concentration, X the latent heat of vaporization, £ water vapor flux 
and q the specific humidity. Overbars indicate an averaging 
operation and primes denote deviations from the mean. 

Raw data were conducted by two dimension coordinate 
rotations [31] and Webb-Pearman-Leuning (VVPL) correction 
[32] to obtain 30-min mean flux data. C0 2 and water vapor fluxes 
could be affected by rain and dew. The abnormal data were 
eliminated following the method used by Falge et al. [33]. In the 
daytime, data gaps were 29%, 8%, 10% and 16% in the growing 
seasons of winter wheat in 2003, 2004, 2005 and 2006, 
respectively. Gap filled data were not used in this paper. 

Light response of NEE 

In most ecosystems, the relationship between daytime NEE and 
PAR can be expressed by a rectangular hyperbolic function 
[33,34]: 



NEE = 



7.P max PAR 
P m!a + aPAR 



+ Rd 



(3) 



where oc is the initial slope of the light response curve (initial light 
use efficiency), P max the maximum photosynthetic capacity, the 
daytime ecosystem respiration rate under dark conditions. The 
model (Eq. (3)) was fitted by the software "Origin 7.0" (Microcal 
Software Inc.). In this study, negative NEE means net C0 2 uptake 
by the cropland and positive NEE indicates net C0 2 emission 
from the cropland. 

To study the environmental factors influencing daytime NEE 
besides PAR, NEE was computed using Eq. (3) and the residual 
(NEE r ) was obtained as measured NEE minus simulated NEE. 
Stepwise multiple linear regression models were used to estimate 
the integrative influences of biophysical factors on NEE r . 
Statistical analyses were performed using the software "SPSS 
13.0" (SPSS Inc.). According to the method used by Carrara et al. 
[17], Powell et al. [18] and Teklemariam et al. [7], air temper- 
ature (T a ), vapor pressure deficit (VPD), soil water content (SWC), 
leaf area index (LAI) and canopy conductance (g c ) were divided 
into many classes and mean NEE r was obtained for each class to 
show the trends in NEE r versus variables. 
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To investigate seasonal patterns of light response parameters, a 
rectangular hyperbola (Eq. (3)) was used to describe the light 
response of NEE every five days and a time series of oc, -P max and 
i?H were obtained. Meanwhile, the biophysical variables in the 
daytime with the 5-day interval were averaged. Stepwise multiple 
linear regression analyses were applied to distinguish the key 
drivers for a, P max and R^. The relationships among light response 
parameters were also investigated at the 5-day scale. 

Light response parameters varied among years. Standard errors 
(«i and e 2 ) were calculated for sunny and cloudy sky conditions and 
total standard error (ej was obtained for both sky conditions using 
the equation as follows: 



e\ +e 2 



(4) 



The difference between sunny and cloudy sky conditions would 
be pronounced if it was larger than the total standard error for two 
sky conditions. 



Canopy conductance 

According to Monteith and Unsworth 
tance (e.) was calculated as follows: 



[35], canopy conduc- 



1 1 



A(R,,-G) + pCpg a VPD 



XEy 



A 

y 



(5) 



where IE is the latent heat flux (W m -2 ), R n net radiation (W 
m -2 ), G soil heat flux (W m -2 ), A the slope of saturation vapor 
pressure to the air temperature curve (kPa K _1 ), p air density (mol 
m i ), C p the specific heat capacity of air (J mol 1 K '), VPD 
vapor pressure deficit (kPa), y the psychrometric constant (kPa 
K -1 ), gc canopy conductance (m s _1 ), and g a aerodynamic 
conductance (m s _1 ) which can be given as follows [36]: 



k 2 u 



(6) 



where ra is aerodynamic resistance (s m— 1), z the height where 
wind velocity measured (m), d zero-plane displacement height (m), 
zO roughness length (m), k the von Karman constant (0.41) and u 
wind velocity (m s— 1). zO and d can be calculated as [37,38]: 



z 0 = 0.13/z 



--0J5h 



(7) 



(8) 



where h is the crop height (m). 



Defining sunny and cloudy sky conditions 

The clearness index is used to describe sky conditions and 
the degree of impact of cloudiness on the solar radiation received 
at the Earth's surface [22]. It can be given by [39]: 



k r - 



S 



(9) 



S e = S sc [l +0.033 cos (360? rf /365)] sin/? (10) 



sm/? = sm</>sin<5 + cos^cosc5cos<» (11) 

where S is global solar radiation (W m z ), S e the extraterrestrial 
irradiance at a plane parallel to the earth surface (W m ), S sc the 
solar constant (1370 W m -2 ), i d the day of year, ft the solar 
elevation angle, (j) the local latitude, S the declination of the sun, 
and oj hour angle. The sky conditions were classified at a half-day 
scale because the days without clouds for the whole daytime were 
rare [39] . In this study, sunny skies were simply identified when 
the half-day mean kt was larger than a threshold in the morning 
and afternoon, respectively. 1.2 times of two-month running 
average k t was used as the threshold so that the number of 
obtained sunny half-days was close to the mean value observed in 
recent decades. Clearness index were plotted against solar 
elevation angles and fitted by cubic polynomials in the sunny 
mornings and afternoons, respectively (Fig. 1). The data falling 
away from the major patterns were excluded [39] . Different from 
the studies in the forests [22,39], ^ observed in the wheat field was 
lower in the afternoon than in the morning (Fig. la). Air pollution 
may reduce the atmospheric transparency, resulting in a small 
incident solar radiation at the land surface. The reduction was 
more evident in the afternoon than in the morning (Fig. lb). 
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Figure 1. Scatter plots and regressions between (a) the 
clearness index (kj and the sine of solar elevation angles 
(sin/), and (b) global solar radiation (5) and sin/ for a winter 
wheat field in April-May 2003. The data were fitted by cubic 
polynomials in the morning (solid line) and afternoon (dashed line), 
respectively. 

doi:10.l371/journal.pone.0089469.g001 
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Results 

NEE Light response and influencing factors 

The relationship between daytime NEE and PAR in a winter 
wheat field is shown in Figure 2. Light response curves in April 
and May were similar but both of them differed from those in 
March. From 2003 to 2006, P max was 46.6±4.0 umol CO, m" 2 
s"\ ot 0.059±0.006 umol umol" 1 and £ d 5.3±0.3 umol C0 2 
m 2 s 1 in April— May, two or three times as high as those in 
March (Table 1). However, the inter-annual coefficients of 
variation (CV) for ot, P max and R d were great in March due to 
large variations of LAI, daytime mean T a , VPD and SWC among 
years. The variations in VPD and SWC resulted from significant 
changes in precipitation in March among years (Tables 1 and 2). 
In March 2005, ot was so small and P max was so large that the light 
response curves were actually close to a line (Table 1 and Fig. 2). 

Table 3 indicates that NEE r declined with the increases in LAI, 
g c , T a and SWC (P<0.001) but increased with the increase in VPD 
(P<0.001). Positive/negative NEE r means measured NEE was 
higher/lower than simulated NEE; or measured C0 2 uptake was 
lower/higher than simulated C0 2 uptake. The impacts of 
biophysical factors on NEE r varied in different months. Factors 
influencing NEE r were sorted as LAI, SWC, g c , T a and VPD in 
March; g c , LAI, VPD and 7" a in April; and LAI, g a VPD and T a in 
May. The effects of SWC on NEE r were significant in March but 
insignificant in April and May (Table 3). 

Shifted trends in NEE r versus variables illustrated the limit in 
linear regression analysis (Fig. 3). With the increase in T a and 
VPD, NEE r turned to increase significandy (P<0.001) when T a 
was more than 25°C or VPD more than 1.1 — 1.3 kPa. With an 
increase in g c , the NEE r decline disappeared, or even turned into 
an increase (P<0.01) when g c exceed 26 mm s~ in April or 
14 mm s in March and May (Fig. 3 and Table 4). No reverse 
trends were found for NEE r versus LAI and SWC. 



Light response parameters: seasonal variation and 
influencing factors 

Light response parameters varied seasonally, with peaks in April 
or May ranged from 64.2 to 86.1 |imol C0 2 m -2 s _1 for P max , 
from 0.070 to 0.087 umol umol" 1 for ot and from 6.5 to 8.5 umol 
C0 2 m" 2 s" 1 for i? d during the year 2003-2006 (Fig. 4). The 
seasonal patterns of P max were similar with LAI but different from 
T a and VPD which have an increasing tendency from March to 
May. Compared with other years, .P max peaked earlier in 2004 
when LAI reached the maximum in mid-April due to fast warming 
in spring (Fig. 4). 

Figure 5 shows the scattered points of biophysical factors (LAI, 
g c , T a , VPD and SWC) versus NEE light response parameters (ot, 
P max and R d ) at 5-day scale during the main growing season of 
winter wheat. Stepwise multiple linear regression models were 
used to assess their relationships (Table 5). P max increased with the 
increase in LAI, g c and T a (P<0.00 1) but reduced with the increase 
in VPD (P<0.001). The factors influencing P max were sorted as 
LAI, g a T a and VPD. ot was proportional to ln(LAI), g c , T a and 
VPD (P<0.001). The impacts of LAI, g c and T a on at were larger 
than that of VPD. R d was proportional to T a (P<0.001). It was 
better to express the relationship between R d and T a using 
exponential equation instead of linear equation. However, the 
influences of SWC on all light response parameters were 
insignificant (Table 5). 

During the growing season of winter wheat, ot increased linearly 
with an increase in R d (P<0.001) (Fig. 6). ot and R d enhanced firstly 
and then declined with the increase in P max . The maximum ot and 
R d appeared when P max was around 50 umol C0 2 m~ 2 s _ . Their 
relationships could be expressed by quadratic polynomials (P< 
0.001) (Fig. 6). 




Figure 2. Response of net ecosystem C0 2 exchange (NEE) to photosynthetically active radiation (PAR) in a winter wheat field, a-c: 
NEE obtained by eddy covariance technique; d-f: light response curves fitted monthly using the rectangular hyperbolic function (Eq. (3)). The values 
of light response parameters were shown in Table 1. 
doi:1 0.1 371 /journal.pone.0089469.g002 
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Table 1. 


Monthly NEE light response parameters (P maX / a and R d ) derived from Eq. (3) for a 


winter wheat field. 










P 

*max 








') z 2 




Month 


Year 


(limol C0 2 n 


l 2 s (umol umol 


(umol C0 2 n 


1 2 s 


n 


March 


2003 


9.2 


0.017 


1.6 




0.388™ 


497 




2004 


15.8 


0.037 


2.4 




0.396™ 


681 




2005 


29.9 


0.008 


0.8 




0.360™ 


701 




2006 


10.6 


0.036 


3.5 




0.431™ 


179 


April 


2003 


40.1 


0.044 


3.8 




0.736™ 


695 




2004 


49.7 


0.067 


5.2 




0.836™ 


746 




2005 


56.4 


0.060 


6.4 




0.806™ 


719 




2006 


44.8 


0.063 


5.0 




0.629™ 


692 


May 


2003 


48.7 


0.050 


5.9 




0.876™ 


402 




2004 


32.8 


0.061 


5.9 




0.619™ 


779 




2005 


64.5 


0.036 


3.1 




0.723™ 


743 




2006 


53.4 


0.070 


6.5 




0.715™ 


754 



P max : maximum photosynthetic capacity; 
a: initial light use efficiency; 

R^. daytime ecosystem respiration under dark conditions. 
Significance of the regression was "™" for P<0.001. 
doi:1 0.1 371 /joumal.pone.0089469.t001 



NEE Light response and NEE-g c relationships under 
various sky conditions 

Light response curves of NEE were determined under sunny 
and cloudy sky conditions, respectively (Fig. 7 and Table 6). The 
analyses were only conducted during the period with large LAI 
(0.8LAI milx <LAKLAI max ) to limit the influences of crop growth 
on NEE light response. Compared with sunny sky conditions, P max 
under cloudy conditions was 19%, 12% and 27% higher in 2003, 
2004 and 2006 but 14% lower in 2005; a under cloudy skies was 
24% and 64% larger in 2003 and 2005, but 5% and 10% less in 



Table 2. Monthly biophysical variables in the winter wheat 
field from March to May in the year 2003-2006. 







T. 


VPD 


swc 


Prec 




Month 


Year 


(°q 


(kPa) 


(m 3 m 3 ) 


(mm) 


LAI 


March 


2003 


8.7 


0.57 


0.184 


42.5 


0.61 




2004 


10.5 


0.69 


0.136 


55.7 


1.47 




2005 


7.9 


0.65 


0.100 


0.1 


1.07 




2006 


10.4 


0.81 


0.106 


0.2 


0.86 


April 


2003 


15.3 


0.73 


0.186 


160.3 


3.20 




2004 


17.1 


0.88 


0.141 


52.4 


5.69 




2005 


17.2 


0.96 


0.139 


27.5 


4.49 




2006 


15.7 


0.79 


0.124 


19.7 


3.81 


May 


2003 


21.9 


0.99 


0.160 


12.4 


2.99 




2004 


21.1 


1.10 


0.134 


46.8 


3.44 




2005 


21.1 


1.07 


0.145 


37.0 


4.61 




2006 


20.9 


0.96 


0.138 


62.4 


4.95 



7 a : daytime mean air temperature; 

VPD: daytime mean vapor pressure deficit; 

SWC: daytime mean soil water content at a depth of 20 cm; 

Prec: total precipitation; 

LAI: mean leaf area index. 

doi:1 0.1 371 /joumal.pone.0089469.t002 



2004 and 2006 (Table 6). On average, P max , a and R d under 
cloudy sky conditions were 9%, 11% and 2% higher than those 
under sunny sky conditions, respectively. However, owing to large 
variation in a, P max and R d among years, their differences between 
two sky conditions were smaller than their total standard errors 
calculated by Eq. (4) (Table 7), indicating unremarkable differ- 
ences in light response parameters between sunny and cloudy sky 
conditions. 

NEE Light response in cloudy days differed from that in sunny 
days. The differences between two light response curves were 
significant in 2003 and 2005 but insignificant in 2004 and 2006. 
Cloudy and sunny light response curves were so close in 2004 and 
2006 that their differences were within the confidence intervals 
due to scattered points of NEE versus PAR (Fig. 7). At the same 
PAR, simulated net C0 2 uptake was 33%, 13%, 23% and 8% 
(averaged 18%) higher under cloudy sky conditions than under 
sunny sky conditions in 2003, 2004, 2005 and 2006, respectively. 
The difference between two sky conditions was more than their 
standard errors (Table 7), suggesting that net CO z uptake under 
cloudy sky conditions was significantly higher than that under 
sunny sky conditions. 

As mentioned above, LAI and g c were key factors affecting 
NEE. After neglecting the effects of LAI on NEE during the period 
with large LAI, the NEE— g c relationships were investigated in 
several radiation classes (Fig. 8). The correlations between NEE 
and g c were described by logarithmic equations under strong, 
moderate and weak radiation, respectively (P<0.001). With the 
increase in ln(g (: ), net C0 2 uptake enhanced more quickly under 
strong radiation than under low radiation. At the same g a net C0 2 
uptake was higher in cloudy days than in sunny days. Nevertheless, 
the differences between two regression lines were insignificant for 
all PAR classes because they were within the wide confidence 
intervals owing to scattered points of NEE versus ln(g c ) (Fig. 8). 
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Table 3. Influences of biophysical factors (T a , VPD, SWC, g c and LAI) on NEE residual in a winter wheat field from March to May in 
the year 2003-2006, estimated by the stepwise multiple linear regression models. 





Month 


Independent variables 


F 




Model 




March 


LAI 


644.33*"" 


0.289 


NEE,= - 


2.854LAI+2.813 




LAI, SWC 


382.77""" 


0.326 


NEE r = - 


2.91 8LAI- 1 7.883SWC+5.257 




LAI, SWC, g c 


284.79*"* 


0.350 


NEE r = - 


2.732LAI - 1 6.507SWC - 0. 1 26g c +5.498 




LAI, SWC, g c , T a 


215.08*** 


0.352 


NEE r = - 


-2.61 1 LAI 1 5.341 SWC-0.1 27g c -0.0297 a +5.496 




LAI, SWC, g 0 r„ VPD 


182.50*** 


0.366 


NEE r = - 


-2.461 LAI-10.812SWC-0.112g c -0.1497 a +1.883VPD+4.499 


April 


9c 


386.84*** 


0.125 


NEE r = - 


0.289g c +2.855 




g c , LAI 


315.42*"* 


0.188 


NEE r = - 


0.274g c -1.004LAI+7.054 




g c , LAI, VPD 


267.37*** 


0.228 


NEE r = - 


0.247g c - 1 .1 30LAI+2.21 4VPD+5.329 




g c , LAI, VPD, T a 


211.68*** 


0.238 


NEE r =- 


0.254g c - 0.97 1 LAI+3.388VPD - 0. 1 55r a +6.275 


May 


LAI 


570.78*** 


0.180 


NEE r = - 


1 .989LAI+8.392 




LAI, g c 


403.46*"* 


0.237 


NEE r =- 


- 1 .738LAI-0.268g c +9.709 




LAI, g c , VPD 


277.50*** 


0.243 


NEE r = - 


- 1 .670LAI-0.264g c +0.752VPD+8.539 




LAI, g 0 VPD, T a 


214.56**" 


0.249 


NEE r = - 


- 1 .777LAI -0.264g c +1 .788VPD - 0. 1 88r a +1 1 .925 



NEE r : NEE residual, fimol C0 2 m -2 s _1 (the dependent variable); 
g c : canopy conductance, mm s~\ 

The meanings and units of T a , VPD, SWC and LAI were the same as Table 2. 
Significance of the regression was "***" for P<0.001. 
doi:1 0.1 371 /joumal.pone.0089469.t003 




Figure 3. Effects of biophysical factors on NEE residual (NEE r ) in a winter wheat field in March (a-e), April (f-j) and May (k-o). NEE 

residual was measured NEE minus simulated NEE. The simulated NEE was obtained by Eq. (3). Biophysical factors include canopy conductance (g c ), 
Leaf area index (LAI), daytime air temperature (7~ a ), vapor pressure deficit (VPD) and soil water content at the depth of 20 cm (SWC). 
doi:1 0.1 371 /journal.pone.0089469.g003 
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Table 4. Correlation coefficients between NEE residual and 
influencing variables when the variables were less or more 
than the thresholds. 





Variable 


Month 


Thresholc 


Correlation coefficient 








Variable< 
Threshold 


Variable> 
Threshold 




May 


25°C 


-0.006 


0.478™ 


VPD 


April 


1.1 kPa 


0.012 


0.268™ 




May 


1.3 kPa 


-0.133 


0.481™ 


9c 


March 


14 mm s _1 


-0.353™ 


0.402" 




April 


26 mm s _1 


-0.375™ 


0.029 




May 


14 mm s _1 


-0.283™ 


0.076 



Significances of the regression were " " for P<0.01 and " " for P<0.001 . 
doi:1 0.1 371 /joumal.pone.0089469.t004 



Discussion 

Factors influencing NEE light response 

In boreal and temperate forests, NEE r generally declines with 
increasing T a and VPD. These trends turn to rise when T a exceeds 
20-25°C [7,17] and VPD exceeds 1.1-1.3 kPa [7,18]. In this 
study, similar phenomena were observed beyond the same 
thresholds in a temperate cropland (Fig. 3 and Table 4). However, 
in the tropical forest, NEE r always increases with the increase in 
7a and VPD without reverse trends [40] (Loescher at al., 2003). It 
may be ascribed to high air temperature and humidity with a 
narrow range in the tropic region. 

Net C0 2 uptake decreased with the increase in T & and VPD 
(Fig. 3), indicating that temperature and water stress occurred. 
Under high temperature (T a >25°C), photosynthesis was prohib- 
ited and soil and plant respirations were great, resulting in a lower 
net C0 2 uptake than the one simulated by Eq. (3). On the other 
hand, VPD controls photosynthetic rate through influencing 
stomatal closure. Under higher VPD, stomatal opening is smaller, 
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Figure 4. Seasonal variation of initial light use efficiency (a), maximum photosynthetic capacity (/'max)' daytime ecosystem 
respiration under dark conditions (R d ), air temperature ( 7" a ), vapor pressure deficit (VPD) and leaf area index (LAI) in a winter wheat 
field during the main growing season, or, P max and Rd were obtained at 5-day intervals using a regression model (Eq. (3)). 
doi:1 0.1 371 /journal.pone.0089469.g004 
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Figure 5. Effects of biophysical factors {g c , LAI, 7" a , VPD and SWC) on light response parameters (a, P max and R d ) at 5-day intervals 
during the main growing season of winter wheat. The meanings of abbreviates were the same as those in Figures 3 and 4. 
doi:1 0.1 371 /journal.pone.0089469.g005 



leading to the decrease of photosynthetic rate [41]. Moreover, the 
decline of photosynthetic rate was likely due to low leaf water 
potential caused by high transpiration rates [42]. Ecosystem 
respiration is large under higher VPD because of high air and soil 
temperatures [41]. Hence, net carbon uptake by the cropland 
dropped. In this study, lower VPD threshold in April than in May 
(Table 4) implied that wheat plants were more sensitive to the 
drought at the fast growing stage than at the senescence stage. 



Around the critical value, the reverse trends in NEE r versus g c 
(Fig. 3 and Table 4) suggested a shift from stomatal limitation to 
non-stomatal limitation on daytime net CO2 exchange. At dawn, 
dusk or in cloudy days, photosynthesis was inhibited by weak 
radiation even though g c was large. Higher g c threshold in April 
than in March and May (Table 4) indicated that net C0 2 uptake 
was more sensitive to stomatal behavior at the rapid growth stage 
than at the slow growth senescence stage. Therefore, NEE r was 



Table 5. Effects of biophysical drivers (LAI, g c , T 3 , VPD and SWC) on light response parameters (P max , a and R d ) at the 5-day scale 
estimated by the stepwise multiple linear regression models during the growing seasons of 2003-2006. 





Dependent variables 


Independent variables 


F 




Model 




p 

r max 


LAI 


267.65*" 


0.812 


Pmax = 1 1 -81 84LAI+1 .061 5 






LAI, g c 


145.05™ 


0.826 


P max = 1 0.3964LAI+1 .291 1 g c -4.6225 






LAI, g a 7 a 


1 1 1 .03™ 


0.847 


Pmax = 8.9964LAI+1 .6064g c +0.71 71 T„- 


-14.4703 




LAI, g c , T„, VPD 


105.74™ 


0.878 


P ma x = 8.1 875LAI+1 .0694g c +1 .96277 a - 


-26.0434VPD-4.6112 


a 


T a 


1 1 .40" 


0.155 


a = 0.001 5T a +0.0273 






T„ g c 


8.62™ 


0.220 


a = 0.001 3T a +0.001 9g c +0.01 51 






r„ g c , VPD 


8.19™ 


0.291 


a = -0.00027 a +0.0030g c +0.0354VPD- 


0.0001 




gc, vpd 


1 2.43™ 


0.290 


B = 0.0029g c +0.0320VPD+0.0005 




a 


In(LAI) 


14.07™ 


0.185 


a = 0.029ln(LAI)+0.040 






In(LAI), VPD 


11.71™ 


0.277 


a = 0.026ln(LAI)+0.023VPD+0.021 




Rd 


7" a 


31.98™ 


0.340 


R d = 0.221 T a +0.897 




Ln(fi d ) 


T a 


32.35™ 


0.343 


ln(R d ) = 0.0277 a +0.145 





The meanings and units of biophysical factors and light response parameters were the same as Tables 1, 2 and 3. 
Significances of the regression were " " for P<0.01 and " " for P<0.001 . 
doi:1 0.1 371 /joumal.pone.0089469.t005 
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Figure 6. Relationships among initial light use efficiency (a), maximum photosynthetic capacity (P max ), and daytime ecosystem 
respiration under dark conditions (R d ) at 5-day intervals during the main growing season of winter wheat. 

doi:10.1371/journal.pone.0089469.g006 
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Figure 7. Response of NEE to PAR under sunny and cloudy sky conditions in a winter wheat field during the period of large leaf 
area index (0.8LAI max <LAKLAI max ). Light response curves (solid curves) were fitted using the rectangular hyperbolic function (Eq. (3)). The 
values of light response parameters were shown in Table 3. Between the upper and lower dotted curves were 95% confidence intervals. 
doi:1 0.1 371 /journal. pone.0089469.g007 
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Table 6. NEE light response parameters (P max , a and R d ) derived from Eq.(3) under sunny and cloudy sky conditions in a winter 
wheat field during the period of large leaf area index (0.8LAI max <LAKLAI max ). 







Sky 




a 








Year 


condition 


(finnol C0 2 m 2 s 1 


(umol nmol 1 ) 


(umol C0 2 m 2 s 1 


) l» 


n 


2003 


Sunny 


56.0 


0.039 


5.3 


0.812™ 


226 




Cloudy 


66.8 


0.049 


5.4 


0.920*" 


302 


2004 


Sunny 


57.5 


0.056 


5.0 


0.839™ 


143 




Cloudy 


64.7 


0.054 


3.2 


0.790™ 


261 


2005 


Sunny 


83.0 


0.034 


3.8 


0.799™ 


231 




Cloudy 


71.2 


0.056 


6.1 


0.881™ 


567 


2006 


Sunny 


57.2 


0.066 


4.8 


0.747™ 


231 




Cloudy 


72.8 


0.059 


4.5 


0.766™ 


325 



Significance of the regression was "™" for P<0.001. 
doi:1 0.1 371 /joumal.pone.0089469.t006 



mainly controlled by g. instead of LAI in April. However, when g c 
was limited in colder March or drier May, LAI became the most 
important factor affecting NEE r in these months. 

The significant impacts of SWC on NEE r in March (Table 3) 
may be ascribed to large variation in soil moisture before and after 
irrigation at the turning green stage of winter wheat (Fig. 4). 
Nevertheless, the sufficient irrigation concealed these effects in the 
following two months (Table 3). 

Factors influencing light response parameters 

Photosynthetic light response parameters were not constant but 
varied seasonally [9,11,19]. For wheat canopy, mean P max in 
April— May observed in this study (49 umol C0 2 m~ 2 s -1 ) was 
lower than the results (62 — 67 umol CO z m 2 s ') obtained in a 
winter wheat field in central Germany [5] and a spring wheat field 
in Manitoba, Canada [43] . Mean a in April— May obtained in this 
study (0.056 umol umol - ) was lower than the value (0.063 |jmol 
Mrnol - ) observed in a winter wheat field in central Germany [5], 
but higher than that (0.036 umol umol -1 ) obtained in a spring 
wheat field in Manitoba, Canada [43]. The magnitudes of light 
response parameters varied among different studies might be 
attributed to the discrepancy in temperature, humidity, canopy 
structure or/ and other factors. 



■Pmax had positive correlations with LAI [9,16,18,19], air 
temperature [17,18] and SWC [9,18] but negative correlation 
with VPD [9,44] . A similar phenomenon was found in this study 
(Fig. 5 and Table 5). Among all factors, LAI was regarded as a key 
factor controlling P max [9,16,18]. .P max increased linearly [18] or 
nonlinearly [16,19] with increasing LAI. The nonlinear relation- 
ship was expected as the leaves shade each other in the ecosystem 
with higher vegetation density [16]. In the cropland, the linear 
correlation between P max and LAI (Fig. 5 and Table 5) illustrated a 
suitable density range for canopy light intercept. 

LAI determines photosynthetic area while g c controls the 
photosynthetic intensity. Owing to the strong link between g c 
and photosynthetic rate, the influences of environmental factors 
(7a, VPD and SWC) on photosynthesis may be ascribed to their 
effects on g c [45]. Hence, factors influencing P max was sorted by 
LAI, g c , T a and VPD (Table 5). Our results in an irrigated wheat 
field differed from those obtained by Zhang et al. [9] who found 
that P max was only affected by LAI and VPD in a dry wheat field. 

In this study, oc was affected by LAI, T m g c and VPD in an 
irrigated wheat field (Table 5), differing from the result reported by 
Zhang et al. [9] who found that a was only influenced by LAI in a 
dry wheat field, a represents weak light use efficiency by the plants. 
At dawn and dusk, the weak light is mainly composed of diffuse 
radiation. Under low LAI, less diffuse radiation was obtained by 



Table 7. Comparisons of light response parameters and simulated NEE (NEE r ) under sunny and cloudy sky conditions. 





Items 


Sky 

conditions 


Average 


Standard 
error 


Difference 
(Cloudy-Sunny) 


Total Ratio 
Standard error (Cloudy/Sunny) 


P max (umol C0 2 n-T 2 s" 1 ) 


Sunny 


63.42 


6.53 


5.45 


5.95 1.09 




Cloudy 


68.87 


1.88 






a (umol (imor 1 ) 


Sunny 


0.0489 


0.0074 


0.0055 


0.0068 1.11 




Cloudy 


0.0544 


0.0022 






ft d (umol C0 2 n-T 2 s _1 ) 


Sunny 


4.75 


0.32 


0.08 


0.67 1.02 




Cloudy 


4.83 


0.62 






NEE S (iimol C0 2 rrT 2 s _1 ) 


Sunny 


-16.02 


1.12 


-2.90* 


1.20 1.18 




Cloudy 


-18.92 


0.58 






For each year, NEE r was calculated at the same PAR using Eq. (3) and the light response parameters in Table 6. The mean values were obtained for two sky conditions 

and total standard error was computed using Eq. (4). 

The meanings of P ma x, ol and /?d were the same as Tables 1. 

Significance of the difference was " " for P<0.05 if the absolute difference between two sky conditions was greater than the total standard error. 



doi:1 0.1 371 /journal.pone.0089469.t007 
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Figure 8. Linear relationship between NEE and logarithmic canopy conductance (g c ) under sunny and cloudy sky conditions (solid 
lines) in a winter wheat field during the period of large LAI (0.8LAI max <LAKLAI max ). Between the upper and lower dotted lines were 95% 
confidence intervals. 
doi:1 0.1 371 /journal.pone.0089469.g008 



the short canopy. When LAI enlarged, the canopy structure was 
beneficial to receiving scatter light from each direction. Less 
diffuse radiation was intercepted by the sparse canopy when LAI 
decreased with leaf senescence. Moreover, at dawn and dusk, low 
radiation led to small g c . Under cold and moist conditions, ot was 
more sensitive to T a and g c than to VPD. 

Generally, the key factor influencing respiration was tempera- 
ture (Table 5). It likely changed to moisture under dry conditions. 
For instance, Zhang et al. [9] obtained that had positive 
correlations with LAI and soil moisture in a dry wheat field. In this 
study, Ra was correlated with T a significandy but not related to 
VPD and SWC (Table 5 and Fig. 5) due to sufficient irrigation in 
the winter wheat field. Furthermore, the effect of LAI on i?j was 
unremarkable (Table 5 and Fig. 5) because R d was composed of 
soil and plant respiration and only the later part was related to LAI 
[16]. 

After investigating the flux data observed in the grasslands and 
croplands over the world, Gilmanov et al. [19] pointed out that 
light response parameters were correlated with each other. 
Because P milx was proportional to LAI, a and R^ correlated to 



P max was the same as correlated to LAI (Table 5, Figs. 5 and 6). 
Positive correlation between a and R d may result from a relatively 
stable critical PAR (PAR under zero NEE) of 110±6, 99±6, 



102±7 and 101 
respectively. 



:9 umol m " s 'in 2003, 2004, 2005 and 2006, 



NEE Light response and NEE-g c relations under various 
sky conditions 

P max in a wheat field was greater under cloudy skies than under 
sunny skies in most of years (Table 6), in agreement with the 
studies in the temperate forests [8,21,22]. Different from the 
temperate area, Zhang et al. [22] observed a reverse phenomenon 
in a subtropical forest: net C0 2 uptake under cloudy skies was less 
than that under sunny skies. It may be owing to much lower 
radiation intensity in cloudy days compared with sunny days in the 
subtropical forest. In the winter wheat field, differences of P mim 
and ot between two sky conditions were insignificant (Table 6) due 
to large variations in P max and a among years. The values of a in 
this study were consistent with those reported by Dengel and 
Grace [8] and Zhang et al. [22] for forests, but different from the 
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results observed by Hollinger et al. [10], Rocha et al. [21] and 
Suyker et al. [24] who obtained a larger a under cloudy conditions 
in the forests and an irrigated maize cropland. 

In spite of no significant differences in light response parameters 
between sky conditions, simulated net CO z uptake under cloudy 
sky conditions was greater than that under sunny sky conditions at 
the same PAR (Fig. 7). A similar phenomenon was observed in the 
temperate forests [8,21] and the irrigated cropland [24]. More 
CO2 uptake by the canopy under cloudy skies than that under 
sunny skies was owing to many reasons. Firstly, the proportion of 
diffuse radiation increases under cloudy sky conditions and more 
light can reach below leaves of the canopies [46] . Photosynthetic 
rates of shaded leaves are promoted by the delivery of diffuse 
radiation [47]. In addition, compared with shaded leaves, the 
phenomenon of saturating photosynthesis easily happens for sunlit 
leaves because shaded leaves often illuminate brightly [48]. 
Secondly, canopy conductance was usually higher under cloudy 
sky conditions than that under sunny sky conditions. Diffuse 
radiation enhances canopy stomatal conductance mainly due to 
the reduction in VPD and blue light enrichment within the canopy 
during cloudy and overcast weather [49,50-52]. Both low VPD 
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